Abstract 


This registration report details the design, 
construction and software support for a digital torsion 
meter capable of recording transient twist of any shaft. 
Graphical results are presented for the delta start 
transient for an induction motor driven electromechanical 
rig and compared to results generated by a lumped 
parameter model for the system. Indications of other’ 
tests that will be carried out and improvements to the 
software and hardware are given. It is projected that the 
accuracy of this device may be improved to 0.02 degrees 
and the signal reveal a torsional signature of the states 
of operation and health of a machine. 


ie 
e 


The Measurement of Twist of the Shaft 
Of gh E1ectromecnanicait Rid. 


A report of the preliminary work carried out in a project to condition 
monitor a crankshaft of an internal conbustion engine 


Author:- John T. Lenaghan 


Supervisor:- G. E. Goodwin 


School of Science and Engineering, 
Humberside College of Higher Education, 
Hul J 

1988 


P ref ace 


This report is. written to support the backdating of 
registration with the CNAA. To this end, it contains 
details of the work completed to July 1988 and some 
results optained since, A httie Ndication iS given of 
the ultimate direction of .the work. 


Jonn 7. Lenaghan 

School of Engineering and Science 
Humberside College of Higher Education 
4tH October 1988 


CNAA Registration report 


Contents. 
Acknowledgements 
Symbols 
Chapter 14 {tntroduction 
-41 introductory 
Review 


4 

1.2 

1.3 Overall aims 

1.4 Application to other systems 


Chapter 2 The torsion meter, principle of operation and 


execution 


Principle of operation 
Mathematical Fheory of Operation 
The construction of the discs 
Fiexibility in design 
The signal generated and the errors associated 
Driving the optoswitches and signal conditioning 
The microcomputer interface Mark 1, ’User Port?’ 
Problems and systematic errors 

2.8.1 Gver Twise 

e,8.e timer Koll (ver 

2.8.3 Variable Interrupt Latency 
2.9 Limitations of the system (BBC microcomputer) 
2.9.1 Frequency of measurements 


2.9,2 Accuracy of time measurement 
2.10 improvements - Mark 11! and beyond 


2.16.4 -Teeneter 760 8 se Dit RISC microcomputer 
2.10.2 Mark il interface 

2.11 Data Processing, Daps 
getbits {recovery of the data} 
daterror fcompensation for faults} 
process fcalculation of the velocity and twist} 
ratios {scaling factors for the coding strip} 


NM MOM PM A PO Pe fe 
OnNAnaewh = 


CHAPTER 3 The Development Rig, Nature and Model 


3.1 The Development Rig 
3.2 The components of the development rig 
3.3 The Mode} 
3.5.4 (he ineguction motor 
3.3.2 The transmission chain 
3.3.2.4. Motor to Drive: Snatt 


jt1300988 { 


CNAA Registration report 


Chapter 3 continued... 


ae. e Drive Shaft to Pulley 
Se 


2. 
3 2,3 The belt transmission 


speed joss 

Torque floss 

Giobal siip 

Model Equations 
3.3.3 DC: generator 


4 The rig mode! coefficients and parameters 
-5 The Mode! Equations tn Summary 

6 Computer solution of the model 

f The computer control of the rig 

3.7.4 Inree Prase Suppiy Controi 

3./.2 Ge@nerator toad 

3.7.3 Software Contro!} 


Chapter 4 Outline of experimentation 


4.1 The start up transient 
4.2 Other tests 


Chapter 5 Fhe Results, Graphs from the Meter and Mode! 


5.1 Comparison between the experimental and theoretical 
5.2 Errors demonstrated by the time domain graphs 


5.2.41 Faitnhniess recording of tow speed ‘transients 
5.2.2 Zones of doubt and sensitivity 
5.2.3 Variable interrupt servicing 


5.3 Frequency domain graphs 
Chapter 6 Concluding discussion 
Appendix 1 Selected Routines in Pascal Code 


Figures, Circuits and Graphs 


Figure 141, Code disc location 

Figure 2, Code disc construction 
Figure 3, Light signa} 

Figure 4, Signal Compensation 
Figure 5, interrupt Servicing 
Figure 6, Flow of information 
Figure 7, Over twist of the encoder 
Figure 8, Photograph of rig 

Figure 9, Schematic of rig 


Figure 10, Belt and pulleys 

Figure 441, Star/detta connection of motor 
General fiow diagram of acquis 

Flow of information 


jt 1300988 | 2 


CNAA Registration report 


Circuit 1 - Opto-switch to TTL interface 
Circuit 2 = thyristor firing circHeit 

Circuit 3 - Firing circuit computer interface 
Circuit 4 - GTO controlled generator loading 
Circuit 5 - GTO computer interface 


Graph 4: Velocity - raw data 


No error detection and the minimum of calculation to 
reveal the velocity 
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The routine for generating the scaling factors is 
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Subsequent processing 
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Symbols 
a angular acceleration 
Cc compliance of the belt transmission 
coefficient of the frictional torque acting on 
rotor 
D d/dt 
e. Young’s modulus for the belt 


f{t} function relating angular displacement to the 
after start 


- function for the rt state variable of lumped 


parameter model 


= tension in the belt arms 
Suffices 
t taut arm 
S Slack arm 


g{i8} the inverse function of ff{T} 


G rotational inductance of the generator 
intermediate results used for Runge-Kutta algorithm 
Suffices 
ron {=<r=<4, 1=<n=<number of variables 
H iteration interval of the Runge-Kutta algorithm 


} current in an equivalent circuit of the motor 


generator 
Suffices 
Swe full model of motor 
x4,¥4 two coil model - stator 
K2,¥2 - "* — - rotor 
D,Q D-Q model stator 
d,q - "* — rotor 


fii current vector for the motor 


Suffices 
3ph three phase vector 
2ph two phase vector © 
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moment of inertia 


Suffices 
r motor rotor, coupling and half shaft 
2 coupling and pulley 


generator and pulley 
currents: in generator 
Suffices 
A armature coils 
- field coils 


integer count of the number of transitions 


Scaling factors for the compensation of strips 


suffices 

{ for the velocity measurement 

2 for the torsion measurement 
stiffness 


SH Shaft 
b belt 
mutual inductance between coils of the motor 
Suffices 
u,v O<U,;V< Ty U<>V 
full motor model 
mM D-@ model, mean tevel 


self tnductance of coils 
Suffices 
See full motor mode} 


number of transitions per revolution at 
measurements are taken 


resistance of co:i;i/s 
Suffices 


ae full motor mode} 


resistance of various components 


Suffices 
9 generator tltoad 
R rotor bars 


S Stator coils 
radius of pulleys 


angular displacement 


Suffices 
general] 
2 


MS motor end of shaft 
Sp pulley end of shaft 
G generator 
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t time elapsed 
Suffices 


{ intervat of ereoder. per ioc 

2 time between transitions at shaft ends 
7 torque 

suffices 

B transmitted by belt 

ai:es 4088 of torque across the belts 

£4 frictional! torque on motor 

f2 frictional torque on generator 

G generator toad torque 

M motor torque (airgap) 

PG terque on the generator due to pelt 

pre initial tension of belt 
Vv phase voltages across motor 

Suffices 

{ ge: full mode! 


ee 
angular velocity 


Kr rth state variable of lumped parameter mode} 

Gs mains supply angular frequency 

{L coefficient of friction between belt and pulley 
ip angie of contact between belt and pulley 
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Chapter 1 
Introduction 
+4 introductory 


"The need to measure accurately the torque, 


or torsion moment, in ae shaft, ... has 
exercised the ingenuity of engineers’ for 
many years."! 


Ford and Douglas wrote this to imtroduce their paper 
in 1949, aware of the fifty years of development that had 
preceded. In the forty years since yet more sensitive | and 
accurate meters applying the currently available tecnno- 
logies have been built. This research is not to develop 
yet another torque meter, but to monitor -the twist by a 
method that will allow application to a number of 
rotating structures and generate signais to indicate the 


States of operation and health. 


The research began as a study into the feasibility 
of condition monitoring a crankshaft wile in service, 
it was taken that the feasibility would be proven by the 
development of sucn a technique. Tins technique would 
measure the torsion- induced im. the . shaft due to the 
torque transmitted and be able to pick out transients 


and characteristic twisting vibrations. 


lrord, Hugh and Douglas, Alan; "The measurement of 
torque in shafts"; Engimeering; 27/5/1949; pp 481 - 
464, 8508 - S07 
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There are many uses for such ae technique, however, 
the initial stages are being carried out on ane electro- 
mechanical rig. This iS to verify the developed 


technique against previous results and theory. 


es Review 


Many methods exist. to measure the torque and power 


output of a prime mover. For example, some momtor a 
Standard Section of the transmission systeme. there are 
other systems developed more recently that will monitor 
the torsional ieee . Of a ane by laser 


techniques. 


The former usually require access for attachments” to 
the transmission chain or even a- standard section to be 
inserted, changing the characteristics of the measured 
system. The latter technique can only monitor the 


velocity at one point from which the torsional vibration 


Cingham, A. J.; "Forque measuring systems. Measurement 
Measurement from rotating shafts, associated 
electronics and calibration methods"; industrial 
Sensors seminar; Advanced manufacturing techniques, 
demonstration heid at Regional} electronics centre, 
University of Hull; 23/6/88 ) 

3Eastwood, y. G. and Halliwell, N. A.; "Laser-based 
measurement of torsional vibration”; SPIE Optics in 
Engineering Management 599; 1985; pp 358 - 366 

ibid, "taser Technology for Condition Monitoring”; 


COMADEM 8&8; Eds. Raj Rao and Anthony Hope; Conference 
at Birmingham Polytechnic; 19-21/9/88; pp 154 - 159 
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nas to be deduced, 


Recent work from the Hitachi industrial laboratories 
measure the torque and power output by the mathematica! 
inference from measurement of speed and the theoretical 
Knowledge of the particular prime movers4. This technique 
~ State observer - has the lhmitation that there must be 
considerable a priori erent of the prime mover and 


the system of which it is apart. 


1.3 . Overall aims 

It iS ar that the monitoring technique 
described below can give information about the particular 
system that is as yet not possible to measure, and this 
measurement may be carried out in a way that does not 
affect the characteristics of the rig. Also that the 
observation. of the system is independent of Knowledge of 
the system, hke the signals obtained from an 


accelerometer. 


The first stage of the project, the stage reported 


40onmae, Tsutomu; Matsuda, Toshihiko; Kanno, Minoru; 

Saito, Keiji and Sukegawa, Takashi; "A Microprocessor- 
Based Motor Speed reguiator Using Fast-Response State 
Opserver for Reduction of Torsion Vibration": fEEE 
Trans. industry Applications 1A-23 (5): 
september /October 1987; pp 863 - 871 
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here, is to develop a monitoring technique. This includes 
analysis of the technique’s imitations and methods of 
deducing the state of the system from the méasurement of 


torsion (and any prerequisite measurements). 


Also reported here is the development of the test 
rig mathematical model which was used in the “first stages 


to verify the results of the measurement technique. 


The principle exploited and tested on this mg iS 
that of described by Pratt?. Tne main difference being 
the use of the newer microcomputers to carry out the 
processing - and display. This report deals with the 
development of this on the drive shaft and the necessary 


computer software to support it. 


The ideas that form the basis of the research are, 
in summar y:- 
A version of the angle of twist detector 


that iS Supported by re) microcomputer to 
allow sophisticated analysis oF its measure- 


ments, 

oPratt, = ee George W. and McMuilin, Paul G.; "Opteo- 
electronic apparatus to generate a pulse modulated 
Signal indicative. of... the state of -Healtn-of the 
mechanical state of a system”; United States Patent 
3871215; granted 18/3/1975. : 

Pratt, dr, George W.; "An Opto-Electronic Torquemeter 
for tagimne Controi”": Society of Automotive Engineers 


SAE paper 750070; 2/1976. 
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A dynamic model of the system that faith- 
fully follows torsion transients. 


1.4 Application to other systems 


When torque is transmitted through a rotating shaft, 
then the resulting twist along the shaft may be measured. 
By the nature of the opto-electronics used the meter 
iS particularly non-invasive, that iS it does not 
require alteration of the system measured other than 
Suitable optical tracks to be attached to the trans- 
mission train and access made for the hght beam. This 
system is useful in areas of hmited = space. By Frey it 
iS possible to use the system on any torque transmitting 


shaft. 
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CHAPTER 2 

the torsion meter, principle of operation and execution 
2.4 Principle of operation 

When any shaft 1S transmitting. power, the torque 
twists that shaft. This meter measures this twist by 
determining the ratio between the time for “two marks of 
Known angular aisplacement at one end of the shaft to 
pass a reference) poimt and the time between the first 
mark mentioned above and ae similar mark at the other” end 
Of the Shaft. The first interval iS inversely 
proportional | to the angular velocity. The second interval 
gives a oa of the aes if ees aia are initially” 


comcident. Hence their ratio reveals the twist. 


if this operation iS carried out Srorumntty within 
a revolution and the times are measured accurately then 
the time history of the velocity, power and twist of the 
shaft may be displayed in adetail. The mathematics 


Supporting this method and its accuracy and sensitivity 


are worked out below. 


A note on terminology; the interval which reveals 
the velocity is referread to as the "perioad’, and that 


which reveals the power, the ‘interval. 
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e.e Mathematical Theory of Operation 


The period of the signal, La, 's used to measure 
the average angular velocity over that particular 
interval. 

D@=( 1)2m. 

(t,) Nes 


tz. 44 
The dispiacement between the two code discs is given by 
(895-84) =(ta)em, 
(ts) M¢ 
(2.2) 
The above formulae are necessarily approximations. 


Usually it is . taken thet .-ang@lier “displacement. is a 


function of time: 


6=fft i 
(2. 3) 
and 
O6<sPF '4-4-. 
(2.4) 
But for the purposes . of this analysis it iS more 
convenient to take the inverse function, g Say, 
t=g{9} 
(2.5) 
and 
S/aqglti=e 1Gict: . 
{2.6} 
This is to allow advantage to be taken of the fixed 
angle over which the measurements would be taken. But 
strictly the following analysis is only certain for 
unidirectional rotation. 
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The time for one period - fram one “of the. codmg.. discs is, 


27 (K+1) /Ne 


: t*-da@, 
{ 
eTtkK /Ns 
(e. Ff) 
=  gfen(Kk+1) /net - gizvk/iet. 
(2.5) 
This gives rise to the Taylor expansion, 
t4= g°’ Letk /no } (2T/ns} 
+! /og’’ fetk/ng} (24/nNg)+.. 
(2.9) 


.The first term gives the equation (2.1) with 
De = Fee 
and and error due to the series truncation of the order 


of the series’ second term (if the series is convergent). 


Simitarty for the twist, the time between successive 
leading edges :s given by 
egrK /Ns etK /Ns 
; a n’ dé - g’ de, 
< | 


G 0 
(2..40) 


where h is the angle function of time for the other 
end of the shaft. integrating out and using the following 


expansion 


to= hf{awk/ng} - gletk/ng} + gO} - N{O}, 
(2.14) 


JtI300988 16 


CNAA Registration report 


defining g(0),h(O) as QO, 
Hhlewk /ngj=glewk/noj +g’ (2etkK/ng} (82-94) 
+1 /og’’ [2mtk/ng} (82-04) f+.. 
(2.12) 
From (2,143: 86a t2. 2) 


tos g’ (emkK/ne) (80-64). 
(2-43) 


From which (2-2) follows ara: the error due to the 
approximation iS - of --the order-eFf the Ani re. term. of 
(2,12). The value for gg’. that is used in oe evaluation 
is waerked .out:.. from  sauation. (2:4). This gives a 
cumulative. error. of 


Vrotg?’ (amk/ng}! [(2m/ng) 4+ (82-984) 41172. 
| (2.14) 


For most measurements the twist is less than et/ng then 
the total error has an upper bound of 


Vyroig’’ L2mk/ng}} (2m/ng) =. 


(2, 73) 
This requires an upper bound on the second 
derivative of time by angle to give an error independent 


of the functional relationship between time and angle. 
This derivative may be expressed in terms of angular 


velocity and acceleration. 


if v is the angular velocity and ais the angutar 


acceleration then 


act = -a 
dec y3 


(2.16) 
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2.3 Fhe construction of the discs 

The ateiteen discs used on this rig are of ABS 
plastic cut to a circumference of 480mm and mounted = on 
each end coupling of the drive shaft®&, To the outside 
edge OF each aisc iS fixed a photographic strip of 


equally tong black and white zones. 


Each Zone. 1. t2nm - jong:  enwine a total of twenty 


black/white transitions at equal angular intervals around 


the disc. To oaumprove contrast the photographic film its 
mounted on an aluminium foil tape attached to the disc!, 
2.4 Flexibility in design 


The dimensions of the rig and early techniques 
developed to mark the discs constrained the discs to a 
circumference of 480) mm. 7 wen. to. .ajiow them. to be 
constructed easily. However, there iS ttle constraint 


on the size that they might be in ae final apparatus. 


The use of photographic techniques aliows a shorter 
strip with finer divisions up to the resolution of the 


interface system. 


SSee Figure 1, Code disc location 
See Figure 2, Code disc construction 
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The signal processed by the computer is derived from 
the changes nm reflectivity at the discs’ edges. 
Considering the reflectivity at one point above the edge 
of each of the discs gives a signa!®S which does not 


follow the perfection of the above equations. 


The errors induced into the mathematical analysis by 
the real avetem stem from the imperfect nature of the 
production of the strips. The formulae “assume that the 
transition occur at precisely em7sniac.Taedian - intervals 
about the circumference of the sions This is not the 


case. The place of transition may in practice vary from 


its theoretical location by *t/_0.25mm on the surface. 


However, any _errors in the measurement due to this 
misalignment of the coding strips will be reflected in 
the measurements as a periodic variation over each 
revolution. This periodic Yariatian may be compensated 
for by the use of scaling factors that are tuned to each 
zone length between transitions. This iS acmeved by 
running the discs at a constant speed and at ae constant 
ears _and comparing the times for each transition 


with that calculated from the periog of one revolution. 


Ssee Figure 3, Light signal 
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ror the signal that 1s used +Or the velocity 


measurement ?: 


feat f3 
where N +S the number - of trans; tions Ze one 
revolution and t-rs the period of that. revolution 
So for: Sach: scat ing factor at 
. oe We ee 
(2. +8) 


Similarly for the measurement of twist each period of 
the original signal giving the interval is given by 


to! ©) ta eg) 


f2.49) 
When the appropriate scaling factors are included, 
Kol titolttKottal sky! ty!=T/N 
i2:. 203 
From @quation (2.18) this reduces to the ratio, 
tol tl /to! oka! ye 
fect 


The scaling factors of att t1mes, pear tits ratro to 
each other. It is then defined that 


Kook i toast An 
| (2.22) 


and ail-other. facters sre scaled to thos. this reduces 
the  fluctwuatrons 1m the raw data vefore processing to a 
residual} bevel which 3S iNdistinguishable from other 


sources oF error. 


eee ee) ee ee 


9See Figure 4, Signal Compensation 
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es Driving the optoswitches and signal conditioning 


{ras iS a description OF the Mark 4 interface 


circuit that drives the reflective opto-switches!? inter- 


faces the signals to the general mwmtertace port. of the 
BEC microcomputer’s (including the Archimedes) usually 
termed the ser Pert, A diagram oF the Sir cat iS 
displayed 12 Ce cin: z, The signals From the transistor 


follow the changes an the . reflected ngnt intensity from 
the -surface moving beneath it. he “test. oF - the. circuit 
compensates for as ramped change in the signal with a 
second transistor and conditions it to Tit levels with a 


chain OF Schmitt triggers. 
2.f The microcomputer interface Mark 1, “User Port?’ 


The - information from the mig consists of an 
electrical signal, ine computer has to fury tris signal 


into bytes ot MForMatkion. 


The interface and communications from the rig were 
developed on a BBC microcomputer Master system and 
upgraded to an Acorn Archimedes. However, this first 


System employs the ‘User Port’ of these computers in the 


et eae ae a a en a a 


'ORS Reflective Optoswitch part no 307-913; Data sheet 
Seis 3/3/eT 
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same active way. The algorithm and use of the- supporting 


component oF this first interface are ascussed. 


Some of the araw bac&ks in terms Of the hmtations 
ana errors 7. data resulting from the use oF this 


mterface are outhned. 


A 6522 VIA supports the User PORT. The features that 


are - particularly used by the . conversion -process -are:- 


A}. The generation of interrupt signals to the system 


Cri on active edges on Signals applied to its 
contro! hnes and 
B) Two general purpose timers which may generate 


interrupts after ae preset time interval 


to. trace the signa! through the. computer system. and 
the software developed, let US consider just one 
transition at the opto-switch. A leading edge appears on 


one of the control hnes to the VIA. 


The internal regsters of tims cmp would hWNave been 


set to allow an ainterrupt to be sent to the .GPu, 


The microprocessor detects the mterrupt status 
$ 


leaves. its current program and carmes out an interrupt 
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service routine!) chosen on the basis-of the VIA status. 
rhe  Tovtine .1or.. ciesring the aAterrupt: due. to :a.. léading 
eacge 9 On the --Sar er ne... 8... ehosern. fms. routine --not..-omy 
cl@éars the srvterruct but -2160--resags a. timer and - passes 
this value into memory and resets the timer for the next 
leading edge. thereby the period of the Signa: - of . the 


cantral} ne 1S turned intoa bytes, 


Trus mechansm iS expanded upon to give the data 
used by the processing program. the general flow 
diagrams!¢ for the data acquisition routines, ‘ACAH1S-* 


outhne the approach 12 more aeta,. 


Acguis interprets the signals Into two groups of 
bytes by the process above and stores them to an onter- 
program buffer for use i a processing program, ’Daps’ 
which finally gives the rig velocity and shaft encoders’ 


relative displacement!3. 
2.& Problems and systematic. errors 


There are two main hardware related errors generated 


by the interface as set up. One 38 concerned with the 

'isee Figure. =., interrupt servicing fadapted from Opre, 

Coiin; tnterfacing the BBC Microcomputer; McGraw-Hill; 
1984) 


YeThe flow diagrams may be found towards the end of the 
document 
'3see Figure -6, Fiow. of -imtformation 
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imit on the displacement measured to an eighteen degree 
arc andthe —-otmer -3 due to .-~— tne... method by which the 


timers are read. 


€.8.1 Over _i wWist 
The encoding discs are given a displacement oF 
between four and nine degrees for the two reasons:- 
nc te allow the recording oF twists in the airec- 
tion reverse to. the direction of motion of the shaft, 
HW)... to prevent the interrupt due to the Start of 
the period from -disc:. -A~ comciding with the interrupt 


for the next . leading-edge When - there. is zero or near 


7ero twist. 


if the twist on the shaft is more than the pneth of 
a black and white Zone. for the arc. remaining after tne 
qaisplacement mentioned above) then the leading edge of 
the Si gnat trem Maretae  ) Cre Of to would come after 


the next leading edge for the end of the per:od!4, 


The acgquis:tion routines copes with iS problem by 
waiting to measure an interrupt of Lhe. -correct nature. 


Nat -+s,..9f 1% Ras Beet or oecessec. “and - mterrupt.. to. -. start 


the period and -mterval timers: then. it wilt not ~ use 

another mlerrupt .ffom. (hat source unt." recewes one 
‘4¢igure ~ 7,.. Over twist ef the ~ encoder 
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that signals the end of the interval. The times 
transferred to the processing program refer to ae double 


period, 


2.8.2. Jiper aoe Oser 

An error arises when reading the 16 bit wide = timer 
eqght Dits at a... tme... The. timer 7. not steppeq when it is 
reaqg and it takes a finite time to read even -nalf of the 
ciocK. To Hiustrate that thus Wn certain circumstances 
iS a problem, consider a timer readmg at the -interrupt, 
Wm | hexadecimal:- : 

FDOd, 

the least -Sianitacan.- byte; O22 - 13. read -mto -one. of the 
Cri regqgsters. This takes three mcroseconds, re) which 


time the timer now readsi- 


ae 2 
as it counts down. fhe -most significant byte, FC is now 
read. When the read value of the timer 1s sent to memory 


it. would: -be- recorded -as:— 

FCO2 
orecisely 266 ..-tde@cwma). off —ine true value at cee event. 
it. +8. Not... pessmr for the acquisition program Oo. guard 
against this. Tne time is sent as read to  Daps, . which 


corrects for these errors. 


Ne 
1 
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2.5.5 Variapie interrupt Latency 

A final error 31S due ‘ee me. mecrococdng of the CPU 
used im mcrocomputers. there is a small variation in the 
times for. the servicing due to the time taken to frst 


detect ihe. signa, the CPU -aonly checks —tor: the  iKQ at the 


end of each machine code instruction the delay mn 
detection may be from zero to LWwo mcroseconds., Tims 
change im agelay 1S reftiectead mn the times read. 


2.9 timitations of the system (BBC microcomputer) 
2.9.1 Frequency of measurements 

phe processing time hMmts the frequency Of nter- 
rupts the system -can fends fence iiimits the frequency of... 
the signals. this imitation constrains the number of 
transitions marked on Seen oe, -e- encoding discs. The 
processing times for the service routines in Acquis when 


runmng> on a BBC Master semes microcomputer are given 3in 


the table below. 


+4 +2 CB i CB2 
QO. 3 13.0 ee eee 3-90 i oe 
entry 23.2 eG. 0 26.6 28.0 
service O9.5 eo: 5 AG, 5 gS, 
other ~ - e565 eG 5 
exit at 4 orT..5 OF 5 ai; Oo 
TOTAL 52.0 66.0 TIS -S 72.0 


where OS. -refers to -the time for the operating system to 


pass the control to the service routines. 
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“entry? iS the time to- decide . on the particular service 
routine. 

service’ :S the time for clearing the condition ana 
carrying out some specific action, 

other’ i he case of the Cer ane << routines. 1..the time 
to send the received data to the main memory. 

"exit" the “time token to resume... the mterriupted activity. 


tA, 2 and: Cat Se refer. te. ne —-source of the  mterrupt- in 


¥ the VIA, 


The maximum frequency of interrupts that the 


cCOMmMpUuLer can cone. Wi ts Oyen by the sum -of- the CBI, 


CB2 service times and time ror addition imitial 
aetection, that 1S! 
194.5 microseconds. 


The maximum frequency that may be dealt with is 

10,4 KHZ. 
Acams places the data from the interface  mto the inter- 
program buffer when the CPU is not otherwise occupied by 
mterrupts. The time required for compiete transfer of 
data collected iS 

3se% microseconds. 
The- time that. must be allowed for processing is 

435.25 microseconds. 


This leads to the maximum frequency of nmterrupts for 
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complete transfer of 
2.5 KHZ. 
A combined frequency 2KHZ is the maximum practicable to 


aliow management of the buffers. 


2.9.e Accuracy of time measurement 

With the -BBC. microcomputer. the timers are clock at 
iMnz. There is a lhmt of one clock pulse on the accuracy 
of measurement of the times. This iS a fixed absolute 
error: particular: to the: system,..of 1 .-mcrosecond. As. the 
periods fae interval decrease = length with Increasing 


speedc the retative “error disoa “mMereases. 


Some digital velocimeters employ algorithms to 
reduce the relative error by the extension of the periods 
of measurement. !/° This approach is not open to the 
current approach as it reduces the frequency of 
measurements. The next section detaiis improvements whch 


reduce this absolute error. 


When dealing with a signal from a stmp made of 20 


black /white zones travelling at 1500 rom, these residual 


rT 


{SBonnert, Richara; "Digital Tachometer with Fast 
Dynamic Response Implemented by a Microprocessor”: 
fEEE Trans industry Applications i1A-149 (6); 


November /December 1983; pp 1052 - 14056 
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errors are: for the Master -semes computers variation in 


latency O-4 us 
clock 4 us 
accuracy ioe. 2 (0.16 degrees) 


It should be noted that the tatency and clock periods are 
fixed SO the accuracy wrt} improve for a particular 
system when measuring at a lower velocity and 


correspondingly worsen at greater velocities. 


2.10 Improvements — Mark I! and beyond 


2. 410,4 Transfer to a 32 bit RISC microcomputer 


With the introduction of a en ane BR tit RSC. 
based microcomputer from ACoOrn, the Archimedes, the 
opportunity was taken to transfer the developed software 
to this faster environment. A tabie above detais the 
processing times for the routines of Acquis. On adapting 
certain routines the processing times within software 
have been reduced to one tenth that § hsted, This left 
enough time  toer--Oaes to “operate im. real time “even. though 


it currently operates in interpreted Basic V, 


The 7YUser Port" Supphed with this computer uses | 
timers that are clocked at eeanz. This reduces the 


absolute error in the time measurements to 5OO 


nanoseconds. 
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To wWipstrate the -lmits...of the hardware used: when 
Geahng with a signal from.-a strip made of 20 black/white 
zones travelhng at 1500: rom, these residual errors are: 
for the . Archimedes variation im 

fatency 3.5-22.5 processor cycles 

ee. 4S 

chock 0.5 us 
. accuracy —-8. 37 (0.04 degrees) 
(Compare with the ontervais for the Master given ins the 
previous section) 
2.10.2 Mark i JN errace 

A large portion oF Acquis is devoted to the 
mampulation oF Various registers within the VIA to 
generate the times. it iS proposed that most of these 
algorithms may Successfully be transferred to hardwired 
iogic,; with the times presented to the CPU directly — on 
its) 6©data DUS. Crew. 2 Scnematic diagram. of part of 
the circ. 

the advantages of tms interface are the ehmination 
of errors im reading due to variable interrupt -. latency, 
the reduction in processing time per signal and the | 
reduction of the absolute error in the measurements § {by 
the use of a clock frequency »>10MHZ). 

USINg the Hustration at the eng of the previous 
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section, the error in measurements of twist will be 0.03% 
(O.004 degrees) 
"eS Data Processing, Daps 

There are four routines that deal with the collection 
data, the correction for hardware faults and the Final 
processing to velocity and Snaft torsion, 
These procedures are named (in this report):- 


getbits 
daterror 


frecaovery of the data} 
fcompensation for faults} 


process {calculation of the velocity and twist} 

-ratios iecating Factors for. the coding strip} 
Except for ae every set of data has to be passed 
through these routines mn turn to be processed. They 
mcorporate the algorithms described above. An 


interpretation Pascal 1S given 18) Appendix :: 
getpits 

Ths routine recovers a 6 helt 6 S6htored. by the 
acquisition program, acquis, in the inter-program buffer 


convert them 


to integer times!®©, 


There are cons:derablie 


differences between Master and Archimedes versions 
because oF the differences in the memory mapping and 
routines in Acquis. A Pascal version of this and its 
106A) flow Giagram for the Master version is given with 
the flow diagrams 
JTISOO9SS a4 
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subprocedures are given!f although in practice machine 


code is. used for -speed. 


daterror 


There are a number of systematic errors introduced 


nto the collected data by the present interface that 
this routine iS designed LO aetect and correct. These 
have been detailed 7.) earher Sections. 

: These errors are detected by this .- routine in. the 


Same way. Four point rolling averages of the periods, C4, 


and pntervals, i>, are Maintained by another routine 
process. Against the averages each new set of data 18S 
tested. 


if there 3S an error due to the over twist of the 


shaft during measurement then the permod measured would 


be approximately twice that of immedately before. The 
period iS halved and time lapsed ar e@ corrected and 
processing continues. 


it over twist continues more than one period the 
interval would be much less than the average. If it falls” 
to below one eighth of the: previous - value it %.. modified 
by a@ding the © prewous period. Ihe following Pascal code 


a a a a ee 


'fappendix 4, getbits 
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Hlustrates the principle elements of ths routine, 


if sample<coldsample[{1})}/8. then 
Sample: =sampletoldperiod{[{1}; 
if (period>1.8*AVperiod) and (period<2.2*AVperiod) 
then 
begin 
period: =<period/e; 
runtime:cruntime+period 
end; 
The ranges on the period are to bracket the differences 
nm the times of the periods due to changes in speed. if 
the interval measured, the value of Sample, iS greater 
than one eighth of the previous measured it is taken as a 
true reading not . subject -to errors. Sy... tms method no 
correction iS required when the shaft returns to the 
normal range of measurement. 
if period>AVperiod+150 then period: =<period-256; 
if sample>AVsample+150 then sample: =samplie-256; 

These. last two conditions of daterror detect the 
clock roliovers white reading. When either tne period or 
mterval (sample) are Jtonger than the previous average by 
some large figure then they are reduced appropnriately. 
The hmit iS set at 150 mcroseconds longer than the 
previous time because of the variation in speed and the 


coding qgdisc strips but to allow the detection of the 


majority of these data faults. 
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process 


Tins updates 


rong averages and recalculates them. 


velocity and code 


equations (4.1) and 


derived from the 


time represents 


runtime, and is used to end the iteration 


ratios 

This routine 
factors for the 
Uniike the -above 
the main iteration 


their services. 


Follows equations 


18see Appendix 14; 
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qisc displacement 
(4.2). For graphical 
accumulated periods 


the time within 


sets up the values 


individual zones 


mentioned routines 
joop . oF =~ daps although 


+O calculate the 


—- 


procedure ratios 
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calculated 


(2.18) .(2.21),(2.22)'%. 
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time 
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scaling 
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CHAPTER 3 
The ~Devetopment Rig. Nature ana Model 


3.4 The Development Rig 

The test rig on which the monitoring technique has 
been developed is described. it will be used for |. further 
work  mto tne ‘nature - of icteenks within induction motor 


drive systems. 


The hst of components below is not i imtended to. be 
exhaustive but describe the principle configuration oF 


the test bed during the . stage of the work reported!?, 


3.2 The components of the development rig 


induction: Meter 
Brook Motors Ltd 
3 phase induction motor R1iG-A31910 


speed 1440 rpm. 
Coil! Voltage 400 /440V 
Coif Current 6A 
Frequency 5SOHZ 
Output Shaft Oe i eee 2 Ee 
Number of rotor bars 34 
Number of stator teeth 36 


Transm:ssion Chain 
Bearings Fenners Ltd 
RHP 1020-20 Seif aligning 


Shaft Mitd Steel 2Omm dia. 
Pulleys Fenners Lta 
SPZ type Cat No. 03120243 
Dia. 160mm 
Bush -“No. 2024 
Belt Fenners Ltd 
SPZ:. T8600 680 for 140/140 


ee ee ee 


19see Figure &G, Photograph of the rig and Figure  9Q, a 
schematic diagram of the same. 
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DC Generator | 
Higgs Motors of Witton 
Mode} W1iO033-1 
Originally wired in compound mode 
Now with 
a separately excited field 
3.3 The Model 
The moae} consists of a system of first order 
ordinary adafferential equations for each part of the rig 
which are solved by the Runge-Kutta algorithme? 


computer. Each jart ot the model is separately 


discussed, 


mMmduction motoret 


WW 
WN 
‘oe 
ne 9 
a) 


The Standard presentation of the electrical 
equations for a three phase induction machine involve a 
six by six matrix connecting the voltages and currents in 
the . three phases for  stater- and. reror. The elements) of 


the matrix are aifferential coefficients with time. 


j gies. 5: eee DM DM DM DM 


iv Ps. 4 
vi DM 1 r tsi omt?s om '4 omM!° om !® ry da 
vel=-lpme} OF er BL omc4 DMS? pme®e eIPGe, 
vet toms!  pm3e@— opm. 3 or 24 COMMS) opm 36 eee 
vil fom’? -ge4) ...ont?> ge “~~ 48. pt” 2 Ge 
ye DM? 3 OM?eS pM 3 pmM24 om °F fr ¥BL £473 
6 64 ge 63 G4 65 .- £ 613.1) 
COThomas, Benku; "The Runge-Kutta Methods”; BYTE; 
_ 4/1986; pp 191 — 210. 
fCeluones, Charles V.; The Unified Theory of Electrical 


Machines; Butterworth; 196f 
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(3,4) may readily be deduced “from the equivalent 


Carcut of the mduction motor but in this State it iS 


cumbersome to solve. Jones®! offers transformations that 
wil} reduce the number of equations in the special case 
of balanced supply. it #s ms model that is adapted. for 


the orime mover. 


bi gpnl={Caltianns. ; 
: (3.2) 


where: (Cg) o'8 29e tPenetormmation matrix from: the . three 
phase currents, {' app} to the resolved equivalent 
currents Of -a@ 446. Ohsase — machine,- {ipo}. The 
ir anstOrmaven.. 16. Sau © 66 tne two phase current vector 
has six components of which two are zero one for the 
rotor system and. 2 ©. er for -the - Stator. By analysis 


it 4s found thatthe = macrix is: 


iC 43 Sia * 
A 
O B 
where B = (2/3)17e O O 
={tfBy' 7" (izeyl fe -o 
-(1/6)'/e¢ -(1/2)17e oO 
(3.3) 
With the transformation: 
[Zapp l= {Cal ' fZgpn}: 
(3.4) 
and dropping the zero rows and columns one obtains: 
{Vv } R ¢+DL O MDcosé@ MDsine@ 
e*te-10' TR 4Ot ~-MOsine MD&ese- | «ti 
v¥1} |MDcose -Mbsiné R +DL_ 0 2pn 
yre MDsin@ MDcose@ of oe ey 
y2 e c (3.5) 
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where iv} are the -voltages across the stator and rotor of 
the equivalent two phase machine and M iS and equivalent 


mutual nmductance between windings (or bars). 


The equations are Stilt non-hnear as they contain 
terms... that. Genend “on . the angular. position -of the. rotor 


with respect to the stator ,@?@ which depends on time. 


Boe aoe { O O ) 
B 0 { 0 O 
O Q Sinwt coswt 
G O coswt Sinwt 
(3.6) 
where w 1s the rotor angular velocity and t the time. 
On applying pGar: 
iv } 2h ii O MD 3 
¥¥ {=| OR ee 0 oe 
oF) | Mw MB! 4t b -t 
2 : = 
O MD Mw i%u = -RE+L -D ia 
Fd ae qd es 7) 
the metor air gap torque is given by: 
Tos C1 AS) 1 ee 
(3s, 8} 


Transforming (it with tea eng {C.) to give the d-aq 
currents then (3.8) becomes: 


T= (4/2) 1 ig ee Peal * 10706) 11 1 ie) tiga). 
(3.3) 
From Jones©!, (d/d@)[L,M}] is 


eee 5, Q Q S1310 37188 . 3196 
mio O fe) sing sine sine? 
G G O sine4 sin@ sine® 
S168. sich 2ine Oo fF oe * 8 
sin@ sine? sine o O O 
e Q oO 


sine Sing sine4 
2 2 
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CNAA 
where Os <= S#<%/3 
G4 = 8447/3 
Mm iS the mean value of the measured 
mutual} pNductance. : 
Cs, +0} 
On : SUDST1 Tet ine ate -f3.92: 


Tye ChqaipDti gig 


where Mis the value taken for the mutual inductance. 
Pay 44 
The currents for tms equation are obtained -from the d-q 
médel given 7. equation eo ae ae Suitably manmpulated to 
make the derivatives of current explicit, ready for 
numerical integration. 
{ 
foe “i R= Moo Mi we OMR (i 
eae 2. eS Z D 
“ ¥ $4 SNS@ 1 MR Mi gw 1 «13 : 
piace 2k cS R rad Gy 
(L - =-M&) 1] -My ML w MR - ® oe 
: ae cy { Ss ee Ss i Q 
-Mv } MR =i we -tL.. tL ow -kiOR 
q S { eee oe q 
ts. t2) 
2.352 The transmission chain 
Do.001. - “MOE 28 Drive shaft 
The arive shaft iS assumed to be isotropic and 
homogeneous therefore iS given a characteristic constant 
Scalar torsional stiffness. Tne torque appned to the 
shaft causes an angular aqdisplacement OF ts ends, a 
twist, in reaction. This iS 
Ts = Koy (8sp-8sm). 
(a, 33 
The Frictional) and windage torques for the drive 
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shaft ana the motor rotor are lumped ana measured 
together. A cubic approximation angular velocity is used 


to model tine Gragg: 
3% cs a aa. De een z ra o E S. 
~4..= G4 to4 as +O 2 (DO) tC4 (DQ} 


ine ecuauen-  governmg. -the moter -at.. the motor to 
shaft couphng, -nence the velocity of the -— motors rotor 
and the displacement of this end of the shaft is:- 


D©@gu= (im-loy- Tey) 414. 


re: oe ee Drive sitet) <= 16 Pulley 
similar to the above the equation of motion of the 


pulley. end: of 4he “shaft is given: by 


D°@cgp=(Tgy- Ty)/!l2, 
(3.44) 


ee We Wee Sees belt transmission©®¢ 
Ths requres a afferent approach from the. shaft 
transmission. A Satisfactory has been derived from a 
consideration of three --aspects:-_af the belt's behaviour. 
(i) The fractional speed loss between the pulleys; 
due to the normal shpping of the belt against the 
pulley Surfaces whiie transmitting torque. 
(41) forque loss - the qtference ve the torque 
apphed at the araving pulley and that sensed at 


A ET RN A a ee 


een diagram of the belt and pulley is provided to help 
comprehension. figure -1C. 
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the following. © this iS taken as due Lo the 
TUCO. 4O6S0S ‘at the bers entry -. and @xit.- points 
Lo the pulleys. 

GH): - Te. Gletat = sip hmit ——. The port... at. which. the 
nei “COSS6s —to- 4f ans torque i. proportion. to that 
apphed because 22 Sshides over the ones surface 


of one of ether. pulley. 


The model makes assumptions about the beit and 
pulley system that nave been demonstrated to be 
inaccurate. oy Cowburn¢e?, The mode! :s not an accurate 


portrayal} oF the belts but gives a Simulation of the 
transmission characteristics when embedded Te) the larger 
model. The short periods over which the model is solved 


allow the following approximations. 


Creep theory?4, where the weit 1S -repiaced by a thin 
elastic. member passed over -flat pulleys, is the basis of 
the mode}, The effects that are considered the most 


import to model are detailed below. 


speed ines 
The elastic belt Stretches In proportion to the 
°3cowburn, David; The mechanical performance of 
automotive V-belts; PhD. Thesis, University of 


Bradford, 1984 

“4cerbert, SB. 6,;. “Force and sl rp pehaviour in V-belts": 
Acta Polytechnica Scandinavia Mechanica} Engineering 
series Mo. 67: i972 


JtISGO9&S 4 4 


CNAA Registration report 


tension set up within ft, However, when transmitting 
torque there is a adfference between the. tensions in the 
two..arms of the pelt. This. ditference causes a .torques on 
the pulleys in reaction. Hence the belt arms are in 


héferent states oF extension, 


AS the Se i. iraveis over the driving pulley it 
relaxes against the dynamic $Pic tie | provided by the 
pulley surface. This friction Supports the aifference in 
tensions. and the reaction torque. The change in extension 
causes the belt to have an exit speed lower than the 
pulley’s.- circumference. Tms process is reversed one the 


following pulley. 


in normal operation, the belt does not globally | slip 
against the pulleys’ Surfaces. rere... is at i¢ast one 
paint on each pulley where the belt iS nela by static 
friction and moves at the same speed as the surface. fms 
pomt 1S - @t-— fF near) the -first. place of Conia. tof beit 
and pulley 18 the agrection of the belt velocity. The 
Slack arm of belt is travelling faster than the _ following 
pulley and the taut arm slower than the cee pulley. 
The following pulley has the arive pulley’s belt exit 


speed and. vice versa. oO On OUTpuUt -Ssnhatt Potates ata 


lower velocity ifiah -2he- mput . shart: This ais the speed 
1oSS, 
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Other effects operate on speed toss, the next most 
important being the radial comphance OF the belt 


aHowing 14-10: "fade out of ime "¥".. grocve used -conA this 


Mg, 
From Cewburn’s, 
Foe 
(DOc¢p-DEcp) /D8ep = Q{Fy-F oJ, 
C3..153 
(3.16) 
Pergue.  --__ 19s85 


There .is- some loss of the transmitted torque due to 
ie Shpping at. the entry and exit points where the 
beit is accelerated to the shape of the pulley, not all 
Of . ine DenGging = Memen..- 7 equred at .entry-<is recovered at. 
exit. The balance is torque tost in frictron.. These and 
other effects are estimated to be, by Cowburne?: 

Tap ee 2 Oe £926. S31 HOPE YRIK} 1 (Fe +h IR 

fer DA, 
for oe S&S SUOSt i toute cO.43 fOr Bb. 

The —-toragge 166s -acts to “reduce the torque im . the 
belt seen by the Following pulley only ang does not 
reduce the torque acting on the belt by the ariving 


pulley. 


Siges! 


The belt sips with -no pomt fixed relative to the 


pulley Surface when the ratio.-of the tensions Te) each 
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arm exceeds the himit given by the capstan 
formula:- 


“Fe /Ps =< Ta geee Ge). 


fs. +S) 
This may reduced to: 
(Fe -Foadmax = CTopatanhtip/es $s, 
is. 749) 
where <f « bajlmas 18: “he reetest. << iterence  'n 


tensions tenable between the tLWwo belt’s arms before 


giobal shppage. 


From. equation... 15.59), the difference im tensions May 
be wrytten exphcitly qaependent on the mstantaneous 
angular. speeds of the driving ana following, 


(Ft-Fs)max=lt/a}efO8qgp - DOFrgi/CIDGgpl), 
fs. 290) 


re ee 


Following the analysis given iy Jones, the torque 
supplied by the generator when Joaded iS given by 


To= DAG(NGSGIrlatkia), 
oe 2S 


with a correction or a the remnant fix. 


This torque acts with the torque apphed from the 
gulley on -the inertia for the generator rotor - resuiting 


a8) an angular acceleration, 
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Pc = tee - Te = Tenis. 


where Tes 
Lee equation 


3:4 7pe -tre 


The detailed 


within the 


model. The 


mode} 1S 


methods 


1S the frictional 
(3423. 
model coef ficients 


described. by 


PHHGUC? 1+-Ore MmOLor 


Electrical 


*S 
R 
L 
M 


Number poles 


Mechanical} 


parameters 
420 volts 
se ohms 
O,.,597 henrys 
QO.565 henrys 
2 


parameters 
0.0421 kgme 


2,86x1073 


3,.47x1073 


Shatt and Belt transmission 


DC generator 


Electrical 
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5 
525 radiansN~'m- 


OOo mM 
1.28%107 Nm! 
G.038c kgm 


parameters 
neglected 
O.:14 
Mee 3 
ties 

po tes 4 

G.25 amps 


torque 


testing OF the 


used follow Jones©! 


registration report 
Sec} 
formulated Similarly 


and parameters 


parameters ror use 


references css the 


and Cowburnes, 
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Mechanical parameters 


| 3 0.2239 kgm 
{Ol Fes 

C5 5. 73x410°% 

ce 4 30020" * 

C7 2.66x1075 

¢g -7.84x1074 

LL 1.4 


3.5 The Model Equations in Summary 


where M iS the value taken for the mutual inductance. 

(244) 
PDi <4 fi vot ol Mow ML ow MR Pat 
D 25 2-5 2 R D 
[Di 2s { =, f+ | -M&u a= MR “ML jsf i) 
Ge ee | os Rr 2 Q 
Di of ate -Mo) poe 1} ML @. MR io ee 
q ee q 4 S ee a : ee q 
ba. 3 [-Mv_} |. MR eee 
d d S { ‘2 1 oR a 


Vso = Ks (8@gp-8gm) - 


: . | See 
Lees C4+coDA+c3 (DA) “+c4 (De) >. 
ta.42) 
D“8oy= (CT y-Tey= Tea 7i 4 - 
(3.33) 
D°A@sp=(Toy- Ta)/loa.- 
(3.14) 
DT = 0.01 (140. 331n[10°EI/R4k ,}3 (F¢4Fo)R, 
tor Ct ee. 
TOR fs ..0, 8 -Seeaetitute 0.4 foro. 
(F -F )j ={1/q}* {De - pa b/ (| D8 |). 
{ 3 wae dp fd dp (3.20) 
(Fe-Fslmax = eTporeatanhtiu/sep). 
(3.49) 
BF 8g. =. Ching - Tes Feels ia, 
(322) 
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3.6 Computer solution of the mode} 


Tie equat 1 One if gecti on 32.3 676 salved by 
Runge-Kutta technigue. 


Gy ptHfp(t, As ne 
GayrHtp(tHH/2, 1464, 4/2) Xe401, 2/2. 
Gy ptHt- 2 +t, x4 tG3 = a) tG2 co} 


ana 


s 


Ar {at Ter) cwxpt (Gp 4t+eiGa p+G3 -3+Gg pr) /6, 


~ 


a Fourth order 


St 
(5 KytGo p/e) 
“4 Bn t93,n} 


where f, 18 (he-€quation oF the eifrerentiat. coefficient 


5 L |. ’ : - , a 
Of xX, if terms of The state variabies: and time, ©% 


3.f Re. eomprter— conta. “e+ the rig 


Experiments wilt be carmed out on 


the states Of 


twist that occur normally and due to "Faults? such as 


Start-up, Star-delta Switching, fast bus 


transfer ana 


Various types Qf loading condtions. ror repeatability 


and flexibihty ihe Swiccmng Of the rig is 


thyristors signalled by computer. 


x ee tparee Phase oupply._ Controal 


controltlied by 


The induction motor nas LWwo configurations Of its 


Stator windings ~t6.- the — three phases of the 


Supply. These 


ae: Star" and “Geka: aie — sister: sontrel. 218s able . to 


rt a et ne nee eee ee es oe ne eS NO oR A 


i : 2 ; ; 
COruLLy anotated flow diagrams are displayed in the appendices 
r41SBtiIngs for the Acorn BBC Master Turbo are AVa: sable Im 


BASIC 1V and 'SO-Pascal, a basting 


iT} BASIC V and ARM 


assembler for the Acorn Archimedes is also avaitable 


JUSOCS&S 
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effect. these Wo... configurations =. and. aiiow:. switchmg = from 


one ta the others®, the Seect:o6n a particular 
configuration ana switching transients . reduces [o.. thet... of 


the SOnirol. of 4S: Seis oF Lnyristors -for each. phase. 


Tins COMTEOF: 1S provided by Che firing crcuts ana the 


Signals Sent to the same from tne computer interface 


3.f.2@ Generator toad 
Experiments . highhght the effects se loading 


Variations on the angie of twist. The joad on the 


generator is sét -by microcomputer by switching gate. turn 


of f thyristors — (6tee) that lie ve parallel resistive 


e1rrecu+ 2s roading the generator armature c® , the 


resistances are chosen to draw from @G to 15.75 amperes 


from the generator. i Ste0s of . O.23A. $nhe contro of. the 


GTOs iS Carr rec ot Or oo tnt eritace> firing circuite?, 


a ee SOTt Ware Contiroal 


COniroa of the motor Switching ana loading of the 


generator iS written Tor the Archimedes only which has 


the capability of interacting mn real time white 


eee ee Nee en et ee Se ee 


CSsee Pour i, star/Delta windings of an induction 


motor 
” : . : : : s 
fefsee Circuits ee ihe -thyt reteor fir iis carcust t and 3, 


The computer interface. 


J P o Se : — ; ~, 
“See Circuit. 4 ‘Load the OC generator 
oS ete : : 

e%see Circuit S&S. -intertace 16 “computer. 
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collecting data. All routines are written 7. HIiSG 


system primitives 


ee 

Gi 

c+ 
a 
@ 


assembler code in the form of - ape 
me) a ’relocatable moduie’ 3° to be used in higher bevel 


programs. 


ee a a NE a ee a A TE A NT AN a 


390Acorn Computers, Archimedes Progrmmer’s Guide, 1988 
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Chapter 4 


Outline — of experimentation 


The experiments carmed Gut On the rig qdescmbed 
above. have two. ams: the development of the software and 


mechanism OF the twist detector to retmeve bett 


a8 


r grades 
oF MmFformawan and the analysis at the effects oF the 
Various  dAAOGeS- oF “operation “Gof “ine ric: we terms - of the 


state + - Iwist-oFf the transmission shaft. These two aims 


6 


coincide in the ~ goal of © continuous “Wenitoring “of this 


system from the ~twist- signais. 


4.1 The Start up transsent 
This transient, which generates some of the severe 
torsional} vibrations, iS used to Configure the torsion 


meter. ang the supporting software. 


The mitial results are compared to the expected 


result from the time dadomain model Ot the mig. Any 
agi¢fferences are carefully noted ana the software 


inspected and improved, 


4.2 Other tests 

To discover parameters that distinguish between 
afferent transient and continuous modes of the system 
operation records. wit. be  -taten of .the state-of twist... for 


each. of the cases. these -will -be anaiysed in; both time 


Oo 
at 


JUISCO98& 
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ana frequency domains 
features. 

Comparisons will he 


determine the mature of the 


JtisO09388 
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‘for significant identifying 
made with the models to 
parameters. discovered, 
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Chapter 5 


The Results, Graphs from the Meter and Model 


‘$ 
¢ Below are the results obtainea from the torsion 


meter and the models = ee rig. Lnese are intended . to. be 
a clear example of the meter’s capabilities by a 
comparison OF the Values obtained with the theoretically 
oredcted. 

These results nave peen collected by the Mark 
interface on a BBC Master microcomputer. 

The. final. SeCtions ~GF the chapter. explain. the nature 
of er errors roe are corrected by = remain after the 


current data processing. 


5.1 Comparison between the experimental and theoretical?! 
A visual mspection of the closeness OT Fit 


between the model and the experimental results. 


A number of the electrical and mechanical 
aspects of the rig were monitored to gauge «tne accuracy 
of the Jumped parameter model. These are presented as 
grapns 14 to 2% and represent the unloaded and 


electromagnetically loaded Start up transient. 


3ithe graphs are drawn after the other figures 
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5.2 Errors demonstrated by the time domain graphs 


be. Faithless recording of low speed transients 
phe Frequency resolution OF the encading aises Ss 
arectly propor tional to the speed, At ful} speed each 
encoding aise Samples the veiocity at approximately 
SOGHZ: ms gwes a Maximum ~<tfansient -trequency of .250HzZ 
recorded. However, at low speeds, tOOrpmM, the fastest 
recorded transient frequency. becomes 147HZz. The resonant 


frequency of the mg is about  2enz, as. is seen m the 
graphs. Compare - OFfanns © ~sn0-.1ii,. the: first. three cycles 
OF twist are lost. and the vibration = Tet: cieaer..wnth 


the -encodinig - discs speed up. 


S.t20 - 4Ones OF > au & and sensitivity, 


ee nes one 


there .are -twe Zones at the -top ang bottom .of the 
angular displacement graphs where the system would have 
to cope with the simultaneous arrival of two interrupts. 
Due — tO. the. -14nterrupt ieatency — He computer -—fads. to: record 
accurate data. The wmaccuracy stems from the nature of 
the timers. In the periogd- between the VIA issGing and the 
Gro acknowledging an IR the timer continues to count 
down. Ths problem wit be mscussed. at. jength. mm. Chapter 
# i 7 ths - extension 10 the time recorded - that causes 
the zones because the CPU would be delayed due to the 


immediately revious ae 


On 
W 
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6.2.3. Variable interrupt. servicing 


jms 1S a detay of ~“petween CG and 4 microseconds on 


top of the - te: recorded: -due.- to. the normal. lateney. phe 
calculated error On the displacement record + /- 


O.08 degrees. 
The width OF the Variation Wn the tail OF the 
gdisplacement Graph 10 158 measured at 0.15 degrees. It is 


proposead that this variation iS gue entirely to 


enna ree 


Ce 


onal 


eeomnestcnnnt Tn a Te ee en a HE RN nee RE CRA ASME HSS, b 


tne Variable delay Wy tne Jee = Service. 


NBN Stata ye —_e 
—_ hae Le BTR, La pee PRIM S OE 


5.2 Frequency domain graphs 

Only. a timmted amount of -anatysis of the  =frequency 
domains oF the velocity and twist signals has been 
carried out.°* the fast Fourier algorithm demonstrates 
the hmitation of the system -by the Nyqust frequency to 
approximately: 2S0H2.- “iis const ent tm due to the- use — of 


20 transitions per revolution for measurement. 


Some sigmficant components) and there harmonics are 
apparent... fromthe: gracn. These _have -been identified in 
frequency. but have not been related sources ow the rig- or 


the processing. 


re es ee 


See goranns. i124: and 13 
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b.e.3 Variable ointerrupt servicing 


Ar eee name = earn eae en ne 


ts 38 a> celayv oF “peiween CO and 4. microseconds .on 


top: of the time. -recordeéed: auc. to... the normal. latency. phe 
calculated error on the displacement record + /- 
O.0& deg ees. 


ihe -widith oF Wie Varigtion... the Lert. OF the 
agisplacement Graph 10-18 measured at 0.15 degrees. it is 


proposed that this Variation is gue entirely to 


the variable delay Wn tne. mterrupt Service. 


Pe Frequency domain graphs 

Only a hmted amount of anmalysis of the frequency 
domains OF the velocity and twist signals has been 
carried out. °& fhe fast Fourier algorithm demanstrates 
the hmitation of the system by the Nyquist frequency to 
approximately: 25SOHZ. -“This constramt is due to the- use . of 


20 transitions per revolution for measurement. 


some  sigmificant components and there harmonics are 


& have been iadentitied nm 


ey) 


apparent ...from ihe ..-grapn. The 
frequency but have. not been related sources in the rig. or 


the processimg. 


L 
te ne eee 


32@cee gSranns:-i12¢ and 173 


cn 
SS 
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Chapter 5 


Concluding _—s_ discussion 


Even at this stage When much expermentation and 


yet to be completed, Signals, after 


i 


analysis need 
elementary Processing, | are generated from the twist Of 
the shaft that are sufficient sensitive and accurate for 
the purpose of deducing the state of the machine. There 
are strong maications that the frequency analysis will 
reveal a Signature that GIVES maication Of the torque 


signature and. the mechanical transfer. to the. detector.. 


AS mentioned these conclusions are as yet poorly 
founded and the remainder of the wark will develop. the 
technique and explore the Signals From the mig under 


different operating conditions. 


The. Work -418 expected —%0- pe— continued to . €xpiore the 


apphcation of this technique to emore general and 
complex . Structures, fOF. S2ampe “a  crankshatt. 


The basis .of the thesis for the degree of Doctor. of 
Philosophy wih) be developed from data collected from = 
crankshaft. There Will be an exploration of the 
feasibility to determine the operating condition of an 


internal combustion engme by the signals from Aan 


t 


encoger, a8 . @escrived  Pere,.- attecned to is .erankshef t. 


CF 
er 
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Appendix 4 


Be ie YR ee See ee of eae Dearest - arcu ae eS Seal aa cet ae 
rc Ge = ee See Me Shee | ee 194) 44076 3 37 Fatt a4 ..¢ 


procedure getbits; 
irecover data from the sideways ram for processing} 
constsraccess=66; 


osword ig Ae de et a 
tyo2e word: =. 0). S35555; 
fudgerec =. packed recora 


case boolean of 
faise: taddr: “byte); 
True 2. -e7At 2 Were 3 
end: 


var dummy: integer; 


Function peeK{acaress: worgs: Byte; : : 
femulates PEEK to inspect a memory location 
from Acornsoft guide to 1S0-Pascal on BBC microcomputer } 
Var fudge: fudgerec; | 
begin 
fudge.int:-se@¢dress, 
peek: <fudge.adar®% 
end{peek }t; 


function samplereturn: integer; {turn data into an integer } 
Var databyte, datavalue, jy: integer; 
begin 
datavalue: =O; 
FOr 21 2 36 2-2 
begin 
databyte:=((-peek(datat+ ))-1) 
+ 296 )M0OD & 4100; 
datavalue:={datavalue) *256 
+tdatabyte; 
end; 
Samplereturn: -<dataval ue; 
end{jsamplereturn}; 


procedure bfrts; {get data from sideways ram} 

Var dummy: integer; 
begin 
params{9]:=readsample mod &100,; 
params[1O0]:=(readsampie div &100) mod &100; 
dummy:=codei (osword,sraccess, params); 
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number: <number+1; 
readsample:-readsampie+3; 
end: {bits} 


beginfbody of getbits} 
bits: 
period: -samplereturn/ibo; 
bi tS: 
Sample:=samplereturn/iEs; 
if (Number=-=e2j)and{(sample>period) then 
begin 
per og: =sSampie; 
SELS; 
. Sample: -=Samplereturn/ie6b 
end; 
end; -getbits} 


procedure fratics; 
Var 4,3 5 4ieeeer. 


begin 
for +2 Sea 
Hegin 
fare 220 te: 20 
begin 
Zorre £1,113: =; | 
if ratiob then zonel[:, j}=1; 
end; 
end; 
if ratiob then 
begin 
readsample:=i120c0o; 
repeat 
begin 
‘ort =} t6- -20 
begin 
getbits; 
daterror; 
process; | 
zone [{O,G}:=zone{O,O}-+period; 
Zone iG, bt} s=zone ld, +} +period; 
zone{1i,O}]:=zone[OG,1i}+sample,; 
Fone. +. 82:2 ene Fs. 1 i4+samp! 6; 
eng; 
end 
until readsample>=36000Q0; 
4 =e ee 
begin 


zone {O, ri= Zonet, + iszone{O,G}/ 
F 


2 
zone (?, + }sZonel 1,4} /zone{t,o} /5 
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end; 
readsample: =O 
numberofsamples: <0. 
ends 
eno. FaLios4, 


procedure daterror, 
{compensates for some hardware errors} 
pegin 
if number >2O then 
OeEGiy 
if sample<oldsamplefi1}/S8 then 
Sample: -sample+tolaperiod{ 1]; 
if (period>1.g*AVperiod) and {period<ede,2d*xAVperiod) 
then 
begin 
period: =period/e; 
runtime: =runtime 
end, 
if period>AVper:iod+150 then period: =period-256, 
aif sample>AVsample+150 then sample: =sample-256; 
end; {daterror } 


procedure process; 
fturn the integer times into velocity and twist} 
Var 3: 1nteger: 


begin 
if ratiob - else 
begin 
-Sample:=Samplexzone[i, 
(numberofsampies/et+toffsetjmod2ed +i} 
period: =periodxezonel[Q, 
(numberofsamples/et+oftfset)mod2edod +1} 
end; 
for -t23- e60nte 3 do 
begin 
oldsamples{[{jti}:-=oldsamplies{j]}; 
AYsampie -cAVSampie+oldsamples{[ 3}; 
oldper1oeds{43+1}:-otbdperiods{ ;]; 
AVperiod _sAVper 1. od+oidperiods{ 5} 
end; 


oldsamples{[i]:=sample; 
oldper:0ds {i}: =per: aq; 
AV¥Vsample:=(AVsample+samp!ie) /4; 
AVperiod: = (AVperiod+period) /4; 
if AVprocB then 
begin 

period: =AVperiod; 
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Sample: =AVsample 
end; 
velocity: =o00/per iod/str spe; 
runtime: =<fruntime+per: od; 
twist: =sample/periiod 
end; {process} 
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Figure 2, Code disc construction 
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Figure 3, Liagnt signal 
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Figure 5S, interrupt servicing 


yy, 
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Figure 7, Over twist of encoder 
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Photograph of rig 
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Figure 9, Schematic of rig 
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Figure 10, Belt and pulleys 
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Figure 11, Star/delta connection of motor 
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General flow diagram of acquis 


Greral flow diagram for ACQUIS, the data acquistlon progrom 
Interrupt service routlnes. 


Interrupt condltlon detected 


ser VIA 


at : 
source’? eS 7 


Yes 


Find source of 
Interrupt 


ol m2 .. TR TIMR2 


Clear the Inter- 


refer to the data sheets on the 622 VIA 
—rupt condltton : 


Corry out read perlod read semple update carry update carry 
secondary actl— end send to ond send to byte byte 
~—vitles memory | memory — 

reset timers 


Ria 


return to Interrupted program 
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Fiow of information 


Information transfer through system 
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CNAA 


interface to computer 
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Graph 1: Velocity - raw data 


No error detection and the minimum of calculation to 
reveal the velocity 
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Graph 2: Velocity - no timer rollover 


The routines for the detection and correction are 
applied to the data for graph 1 


1600.00 
1400.00 
1200.00 
1000.00 
800.00 
600.00 
400.00 
200.00 


jt1300988 79 


0.00 


2.00 Time/seconds 


0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 1.380 


0.00 


Veloclty/rpm 


1800.00 


CNAA Registration report 


» 


Graph 3: Velocity - as previous and no over-twist 


The routines for the detection and correction of over 
twist are applied to the data 
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nh 4: Velocity - as previous and encoding disc 
compensation 
The routine for generating the scaling factors is 
applied to the data and those factors are used in the 
stihseatant ornraccinan 
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Graph 5: Velocity - Model 
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Graph 6: Zones of doubt 


The areas of the graphs where interrupt latency 
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Graph 7:Angular Displacement - Raw Data 
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Graph 8:Angular Displacement - no timer rollover 
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Graph 9:Angular Displacement - aS previous and no over twist 


2.00 Tlme 


o) 
rr 
s 
oe 
Oo 
@ 
a 
—" 
| oO 
: ee 
s 
wo 
oO 
NI 
‘a 
a) 
Cm a 
ee 
w= 
xt 
oa = 
= : 
e Som) 
re) 
© 
o 
wet 
ro) 
= re) 
i . 
o 
tom | 
oO 
e 
c 
= oO 
= Te) 
e 
Oo 
a ss 
iss 
e 
ae : Oo 
sama SE 
Fee ee 
Ce x 
e 
cnmemngpinte o 
e 
_ 
~ P 
~ | 
ro) 
i Bee ee ee ee SS ee eee 
° 3 ° 
© Oo r=) oO Oo oO 
et 
Oo oO oO roo) oO ro) oO 
etn = j A — Oo oO = ro 3 
af WN e ® e e 
a = a ook i se rel wo tb “I = 


jt) 300988 | 86 


CNAA Registration report 


» 


Graph i0:Angular Displacement - as previous and disc 
compensation 


\ 


| 


— 
& 
a = 
et ' 
@ 
-“< 
N 
— 
e 
6 
e eee ee ok eee = ; 
o 
: o ao Lew | ow | oO 
~ 
ao Oo (es) i | Oo oO ep) 
ay = es ® ® oO oS oO oO Oo a 
af N 6 ® e @ ® 
= = “ = pave = w t N oO 


jt! 300988 oF 


2.00 Tlme/eeconds 


0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 1.80 


00 


0. 


Displacement /degrees 


138.00 


CNAA 


Graph i1:Angular Displacement - 


16.00 


jtl300988 


14.00 


» 


Registration report 


Mode} 


Lemes esevgiArAs = 


——— ee ee 


0.60 G.80 1.00 1.20 1.46 1.60 1.80 2.00 


0.40 


Gg. 20 


8.00 


} 


0.00 


. 00 
4.00 
2.00 
0.00 


88 


uniform 


TONERS Re TRS a eam PRT 


Linear 
89 


era enn 
HPetebery Moen eUOH te 6a ty tytn tn 4. ta dae yHrty HRS T+ Ls NE? HO Op AEDT On bo REY DY hor ey enEsET TY Rete Pyne ty tOEs Cyto MEbetetet Dt: eoLert; fre BORGER Fave + ry Teg HF HPET TY REE Pvt oeyNy Be te COTE tre eg ey DyRoane fray: f 
weedy tel 


Registration report 


of the 


velocity section 


be 


Frequency analysis - velocity, 


e 
» 


CNAA 
Graph i2 


OOTY Cette roe, teay kr wy Ha Re syere Oy Heit f Ge wr te CORTON DC Dy Rpey MOT rte be enOTtaRtRT trac TmD ya MT eters BH Paty ty ny Orbe bate Rrty Pye Oey Opty RATT tarteee ty Rete Habe oOr tad oe Werebetpteter st omy DubeteerTyerr Tt, teen Dy y aubT by toneBs te OY Ey M9 Hy ROeY Er eres ty ORs tote Pity? Meher (Phy Reve Ee ryt mare; 
Aa 


Peeve Gees te WELT) few wEdt yhe CpepeeRERE eT hye oe tede ty Or) tm 8: 


30eh OR Ty Spans yey temeeT HE bs SroN ey fa ow Ty te: A20GT7 Hew Ute oye= cory Nols NERS tee ATE Tre OPLeT mee eE re trts Rate NEM Rye bs ekety beds Cry Or’ 


UO ACH tty tana hy Dae re rte to PEt Sa Own be year Rear BE Hm aFTY Hy OR OD Myer SOP TEAT Oat Cte Het tyt REDE ereG Th Ort; Wet g Orde Net 7MRapeRtc tetwerte rts Erte tats ze Tw 


WR te mean fe enn eo We 
stryee tem rete een ee ee: Arete ERO ERE AERA A are AMAA sh ge HATS rRN PT sali te HS 


jt} 300988 


| 
j it 
themnretmen erred ee 


{. 


4 
PAU ah a ak eS ET ae Dn ie Gee et omni nt aching al ‘ F 


1. 
in 
Tati 
scrmasta wevtate yee aatetas estate rm Dyn mated ae AEE PTE MPO Fated HS ty tomes orto: QC) 


PRR emery radar he ty treet yt : ican afprsemrnaneccnmmmnet ee sealed Soneiiaeanenmn ‘ena renee t 1? ramemnens cacao ene M 
8 - tia) 
ann] Ww A 


Novomatidoct 


FFT af 197. Samples 


Hl 


ea 


CNAA 


~~ 


Graph 13:Freq. 


Novmaliced 


4 
' Seer Se escamod aenewsnsiffmsancne foumnncsaeesifmes narnia sovesansenlfneseaaasnrnftmamannnnaif eamconesece ffm sein f 
t 4 ie 
Se ecal 
As 69 MAD oe 


at AVAL somles 
oewoge. volul ami 


jt1 300988 


analysis 


Wem Po te ro tn De) 


tg 9835 m 


ee hk Oe 


Pete Cr LO oe ee be Le 2 aah 


eh mb eek tae Meee 


Pe rr 


Ub dares romenny t: ied, Hempendsts mtr, Om oh: had 
‘ 


4h 


Oph, bath ote bo hanes 


ee bee 
fis 


: 
Pica datikin. 


Registration report 


uniform velocity section 


: = a 
A = -Y 2 
2 4 2 FY, 
3 : 2 a 
| 3 ra . 
z 2 z 

= yu 2 
= z Ff g 
H z 3 Ps 
3 a= ig 
: ? bi < 
5 4 4 

aie rd A RS SES ae, 


90 


. ‘i ‘dO SLNIWNYLSNI 3JEVNVLVM 


spucosees/eui, O00°2 08°. .. OS"! Ov “l Uc 4 a) aie 08s-o 03°0 Ov “0 02 °O aes 


oi 


eat: 6 


motor, 


HeaS [= 


Registration report 


No load 


mode } 


acceleration, 


e3s° I 


Verification of mode} 


3a! ¢ 


CNAA 
Graph 14 
jt1300988 


: | cag°Z 


Ud ,P_OLETESIY 


a ; ‘dYHOD SANAWNYLSN! SAVNVLUW cs 


pisowh sees; GO°2 abt bet Get Get. eet. | ee 8 noo oes OES 00-0 | 


94 


ear °Z- 


motor, 


rAS‘{- 


electromagnetically 


Registration report 


bau G- 


acceleration, 


mode } 
Loaded 


cau &* 


Raa} 


Verification of model 
N 
Te} 
me | 
q 
op) 


¢ 
° 


est *2 


CNAA 
Graph 14 
jtil 300988 


S446 


uo,pPoOLeTes5y 


CNAA Registration report 


~ 


Graph 15:Verification. of modei - acceleration, 
experiment 


motor, 


No joad 


uap ¥99F 
ies 


¥ 


i 7 Nag 
BY ww 6 tbe heels as. 


aw * ys fe 
a eee ree es 


nm 


oS on -%> see 
a: . 


4 


ats ood 


here ama 


Los. 


Pha 
tf 


» 


< a ey ¢ 
Bis 


% a 
Dallas) tyabebs 


5 a 
Senne: 
a? 


oe 


A ’ 
inpae pt ne, ae 


URLS Rel oe 
‘ee? atts a 

oa OOS Wee el oe bed 
> a eh oe 
A 


gees 


FRAT 


jt1300988 92 


aft 
a Ad 


tri5¥ by 


Cd 


«aes 


rt 
o.. 


pa c= Sa Stn goer aR st ot 


+e Ops; 
: a este eee 


oa a | * 


Heat 


Nae 


yy 
Ye 


” 


, 


Sy 


wae 
we! 


city 
— 
. 


0 


Mm 


s 


Veloclty/rpm 


1800. 


CNAA Registration report 
Graph 16:Verification of model - velocity, motor, model 
No load 
= = = = . 
Oo oO ~ -) Oo oO a) oe 
de) a N ow} om) | © Lem | 
—“ — — — a a ~ 


jtl 300988 93 


2.00 Tlme/seconds 


0.60 0.80 1.00 1.20 1.40 1.60 1.80 


0.40 


Veloclty/rpm 


1800. 


CNAA 


Graph 16:Verification of modei - velocity, motor, model 


1600. 


jtl300988 


~ 


1400. 


1200; 


1000. 


Registration report 


Loaded electromagnetically 


2.00 Time/seconds 


0.80 1.00 ae 2 t 1.40 ee 3% 1.50 


0.60 


0. 40 


800. 
600. 
400. 
200. 


93 


E] WATANABE INSTRUMENTS CORP. 


Registration report 


HOISTAID 4oreu tad wag 2op aft 


Ne 
ae 
we 
CS 


x ek SRA ha Oa-..#° 
hi tig! 4 


Sat 


Des 


94 


CNAA 


i. 
O 
rw) 
Oo 
& 
4 
+ 
~ VY 
= 
ee wa 
OS, ts 
oe  & 
—- OF 
oo x%..¢ 
> ee 
| te 
, # 
O 
O 
= 
Gn 
@) 
c 
oO 
ae 
G 
OY 
q. 
wv ; 
- i 
hs 
<. 
O 
40) 
t 
© 


jt1300988 


LIOISTAID JOFeW Jad wWdug' 2ov rOfeyS 


i. 

a) - 
— fu 

i o 

T] +7 

i co) 


Ay) 4~ 
~s pee OO 
hs a 
ti) $07 = 
ee sic 
ae 

CO 4 
— OL 

yx 

> Samant 9! 

} 
3 


MAA 
aph 


_—/ 
Gor 


et 


Spu¥csesseowy] ehh eae gs 


stration report 


Regi 


CNAA 


mode } 


motor, 


48: Verification of model 


Grapn 


current, 


Inad 


NO 


O08 


Oot 


Oy i 


US ag 


GU T- 


00° 


00 ct 


WT oo oz 


| 00°82 


O00 °9¢ 


H/}uedins esvyd 10yOY 


95 


jt1300988 


spuosessewt] O02 


report 


ry 


rstratio 


eq 


FE 


CNAA 


motor. 


of model - current. 


,oOn 


O8‘t the me Oy | Ue 


oot 08 °0 090 Ot “0 


‘dHOD SINAIWNYLSN!I 38VNVIVM EJ 


02°0 > 00°0 
“0 *9s- 


bo Ba 


mo ot 


oo: F- 


00°? 


O00 °Cc!l 


00°82 


OG Ss 


H/pueisinsa esvyd 1040p 


? 
} 


> 
} 


$3009 


+t 


— 


od 


LIOISTAID JOPEW 


eo. 


eee (oe ed ee Coke 
4 


= , ' 


re ar 


Pijloa eee 
z SI sey* ej sot: 


acd 


t 

tf 
~. 

z 


= 
ic 


arr 
“Dp 
| r 


- 
~ 


C 


“¢ 2." 


~— 


wet. +2 eee 


Pal; 


+. 


O }- 
, 

4 { 
| 3 x mA) : 2 Scene ae bs 
a es Sea ES ay ioe, tn Aa Sie a 


Ayer yp mee 


Gewese eB are 


} 


ot We ee os 
2 iid 


oe A Pe ey Pe Re 


eo 


Ver} 


‘ 
e 
liad 


“yh: 


49 
ype 


ot eet peers te & 


Graph 


‘ to ~ ie 


YO, ia7eos 


‘at! Oe ae 


(ey i j ‘ ae oh 2 ARS ce 2 vee mee) 
Fe ce Pe Ee ee GA Ot OR wh et a bade SEK meee Oe nO - se ree we ee een ew eed ie are 


ion report 


LIGISIAID JOPew 4od VOL sa[LOS 


lier 


a 
3 


atid 


are v 


ae, 
Hine 


* 


trat 
motor, 
t 


egis 


B * 
Ley 
aw 


ee 


awe ed 


re 
enper imen 


current, 


heh» 
i] 


APT SS hee 
A 


j so ; . f ; Wee EE 
Ai oe: hed Sas ar eee aNd fee, 
J . 


tion of mode 


;Cca 


* 
a) 


Pe ig 
5 


Yer | 


ig 


+4 

Sat $3 
iste te dil 
v@ 


CNAA 
Grapn 


iG # Sar hy 
at tes ra ty tee 


. 
A445 s wf bey 
4 ae 


eenadnaad cet we 


96 


jt1300988 


spucses/eowy] 


Registration report 


CNAA 


generator, 


current, 
mode } 


Graon 20: Verification of model 


00 ..< oe 4 ke Me | 


a a i ee re 


Ov “1 


5 yaaa § 


U0 t 


‘dHOD SINSIWNNYLSNI J8VNVLIVM 


Ov “0 02 0 00°0 
03000 


ede 
DA00*Y 
300° 2 
0300°8 
400° 


Pag. ot 


aus | 


eer okt ae ah 


H/};UeisnNs iozpoiauser 


97 


jt1300988 


generator, 


Registration report 
experiment 


current, 


tion of mode} 


1ca 


] 


Ver 


e 
o 


CNAA 
Graph 21 


imeem  Jofem,. 180 Ve, Iles 


, 


“earn? 


a er 


eS °; ane 
Hyg AF pas 


‘ 
~ 
= th Sao ie 
os pv ere ee 


7 
te 


yar 

2 ¢ 

Ps NEY 5s 
LeAZ 


So it 


gages: 


ae 


= 
_ 


a 


a: 


ari 


os 
SEE: 
3 LA te 
Soe 
48> 


wie 

24%: 
bo 

ye % 


<a 
ee 
x? 


ie 


2 


Nes 


42: 


b ¢. 
2 aoe) 


1 om 
ve 


—- eae 
ar Fo Aro ee cee g a 4 
st ep: Sri 

Y 

AY 

~ . 

AS See pe, 7] 
eas ‘ 


? 
ery; 


Se ee ey | 


et) 


7 eh om 


98 


jti3soogss 


program daps 


ST) aa 


def ault 


set up system, 


assign program 


vartables 


get screen 


display title. . 


screen 


getkey 
walt for a Key- 


—press 


selector 


enter menu 


driven 


routines 


oe 


procedure matn 
(_) 
molnintttal 


Inttloltstng 
procedures 


malniterated a 


tterated 
rout ines 


malnconclude 


concluding 
rout [nes 


EXIT 


rocedure. matninittal 


cepture 
load and run 
acqul silion 
program 


Yes 


store 
recovered 


dot a7 


storage 


store recovered 
rew date 


Yes 


disc data? 
recover 


get data from 
a dise file 


inttlal 


set vorloebloes 
for following 
procedures 


plot Yes 


velocity 


2 vertSC 


plot the scales 
for velocity 


Yes 
vert SC 


No plot the scales 
for code disc 
di sp lacement 
kor 250, detoelat 


plot the 
horizontal scale 
ond dete/t!me 


freph on Yes 
VDU ? 


draw axes and 
scales on 
screen 


EXIT 


grophscreen 


frocedure matniterate 


ENTRY 


Yes 


getbits 


get one set of 
date to pro- 
cess thls go 


‘Yes 


slit gl 


one tteratlon of 


No system mode] 
error Yes 
compenseat | on det error 
7 
No | to correct for 


errors generated 
by the hardware 


process 


convert data to 
velocity and disc 


displacement 


Yes 
data out’? 


ft leout 


False (o> 


Plot 
veloclty 
? 


Yes 
groephp lot 


send veloctty 
to plotter 


Yes 
grephp lot 


send code disc 
di sp lacement 
to plotter 


plot Yes 


on screen 


? 


vdup lot 


send velockty 
ond displacement 
to the screen 


EXIT 


monttor 


monitor the 


veloctty and 
displacement for 
each Iteration 


procedure malnconclude 


ENTRY 
rundown » | 
test? rundowbpert 
eost squares 
ee f lt 


Uses the plotter 


routipes for 
wie the port , 
sending the ‘home’ 


command end closing 
the port. 


command to 
plotter 


Wet f ord Electronics, 
Durpout J&& ROM 


of screen routines 


ose files Filing system 
(*elose) enicn 


procedure getbits 


ENTRY | 


blts 

get next 

set of data 
samp lereturn 

convert to 

pertod 

bits 

get next 

data set 
somp lereturn 


convert to 
somp le 


of semples=2 
and 


. False 


set period 
equel to 
sone le 


get next 
data set 


convert to 
somp le 


